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Dark matter annihilating in our Galaxy's halo and elsewhere in the universe is expected to generate 
a diffuse flux of gamma rays, potentially observable with next generation satellite-based experiments, 
such as GLAST. In this article, we study the signatures of dark matter in the angular distribution of 
this radiation, in particular the deterministic ones. Pertaining to the extragalactic contribution, we 
discuss the effect of the motion of the solar system with respect to the cosmological rest frame, and 
anisotropics due to the structure of our local universe. For the gamma ray flux from dark matter in 
our own Galactic halo, we discuss the effects of the offset position of the solar system, the Compton- 
Getting effect, the asphericity of the Milky Way halo, and the signatures of nearby substructure. 
We explore the prospects for the detection of these features by the GLAST satellite and flnd that, 
if ~10% or more of the diffuse gamma ray background observed by EGRET is the result of dark 
matter annihilations, then GLAST should be sensitive to anisotropics at the sub-percent level. Such 
precision would be sufficient to detect many, if not all, of the signatures discussed in this paper. 
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I. INTRODUCTION 



The spectrum of gamma rays generated in dark matter annihilations has the potential to serve as a 
valuable tool in exploring the physical and astrophysical properties of dark matter. As the dark matter 
annihilation rate scales with the square of the density, one would naively expect the brightest sources 
of gamma rays to be nearby and high-density objects, such as the Galactic Center dwarf spheroidal 
galaxies , or other local dark matter structures Q . 

It is not necessarily the case, however, that these sources are the easiest targets to detect the flux of 
dark matter annihilation radiation. In particular, the distribution of dark matter in the Galactic Center 
is not well known. Although N-body simulations predict that dark matter halos should contain high 
density cusps in their centers 0, Q , mechanisms have been proposed which would modify this prediction. 
This is particularly true in the innermost regions of these halos, where additional phenomena involving 
baryons or central black holes may be relevant @ . As the rate of dark matter annihilations in these halos 
depends critically on the density in their innermost centers, it is extremely difficult to reliably estimate 
the brightness of dark matter annihilation radiation from this class of regions. Although the profiles of 
dwarf spheroidals are better constrained, to a minor extent these targets suffer of similar problems. 

Astrophysical sources of gamma rays can also provide a formidable background for dark matter searches 
from local, high density sources. In particular, the Galactic Center has been found by H.E.S.S. and other 
gamma ray telescopes to contain a bright source of very high energy radiation Q. The presence of this 
source poses a very serious challenge to future dark matter searches in the Galactic Center [1] . 

With these issues in mind, it is difficult to assess the expected signal and corresponding background 
expected from annihilating dark matter in nearby, high density regions, even within the context of 
a well-defined particle physics model. Moreover, even if such a source were observed, one might be 
concerned whether it could be reliably identified as dark matter radiation, rather than as another class 
of astrophysical object. Only with a high precision measurement of the gamma ray spectrum, ideally 
including the detection of mono-energetic lines, could point sources of dark matter annihilation radiation 
be conclusively identified 

An alternative target for indirect dark matter searches is the diffuse, or unresolved, contribution to the 
gamma ray spectrum 

[iS[ni,[il- Cosmolo gical data — in particular the cosmic microwave background 
anisotropics and the observations of the large scale matter power spectrum — strongly favor (relatively) 
cold dark matter. In turn, this implies that structure formation takes place in a hierarchical manner, 
i.e. small objects formed at earlier times and larger objects formed at later times from their subsequent 
merging. As a result of this process, dark matter structures ("halos") are expected to host large numbers 
of clumps ( "sub-halos" ) , which are the relics of the formation history. Given the fact that the annihilation 
signal scales with the square of the dark matter density, the presence of substructure can result in an 
enhancement of several orders of magnitude in the gamma ray flux with respect to the naive estimate 
for a smooth distribution, as first noticed in [Tsj . For a simple argument justifying this conclusion, see 
Appendix [X] 

The gamma ray signal from dark matter halos (including our own Galaxy) depends crucially on the 



2 



fraction of sub-halos that survives the hierarchical merging process, as well as their density profile and 
spatial distribution. The lack of knowledge of the dark matter power spectrum on very small scales, the 
difhculty involved in calculating gravitational clustering in the deep non-linear regime, and uncertainties 
in the effects of the tidal stripping of sub-halos by both dark matter and baryonic structures conspire to 
make current estimates of the average annihilation rate uncertain by orders of magnitude, independently 
of unknowns regarding the particle identity of the dark matter candidate. At present, it is even unclear 
whether the diffuse dark matter signal should be dominated by our Galactic halo or by the extragalactic 
component. Furthermore, even if our Galaxy's halo dominates the diffuse gamma ray spectrum, it is not 
clear from which direction, and to what extent, the emission from dark matter sub-halos will dominate 
the overall flux from dark matter annihilations. 

A lot of attention has been given in the literature to studying the spectrum of gamma rays in order to 
identify signatures of dark matter (for a recent example, see Ref. @). Although this has the potential 
to be a powerful diagnostic tool, it will be challenging to clearly separate a dark matter signal from 
alternative astrophysical sources unless very precise measurements are performed with a larger number 
of events and/or over a large energy range. This is especially true in the case of diffuse radiation and 
in scenarios where dark matter annihilations only provide a subleading contribution to the total gamma 
ray flux. 

In this paper, we mainly focus our attention to the complementary information provided by the an- 
gular patterns in the diffuse radiation generated in dark matter annihilations, stressing in particular on 
deterministic features. In Sec. [Til we provide a simple parametrization of the Galactic and extragalactic 
diffuse dark matter fluxes and discuss under which circumstances each of these components dominate. 
Secondly, we summarize the key signatures and distinctive features of the dark matter diffuse gamma 
ray radiation if dominated by the extragalactic (Sec. IIIip or Galactic contribution (Sec. lIVp . We also 
discuss the prospects for the detection of these signatures by the forthcoming generation of instruments, 
in particular the GLAST satellite (Sec.|V]). We summarize our conclusions in Sec. IVII 



II. GALACTIC VS. EXTRAGALACTIC DOMINANCE 



In this section, we compare the relative importance of the Galactic and extragalactic contributions to 
the diffuse gamma ray emission. Depending on which of these dominates, different signatures are expected 
in future gamma ray observations. Although the average dark matter density in both cases is pretty well 
known, because of the effect of substructures discussed in Sec. I and illustrated in Appendix |3 the 
very prediction of the average flux is a highly non-trivial problem. In general, apart for particle physics 
details, it depends on several astrophysical and cosmological parameters. In the following, we shall 
introduce simple phenomenological models which condense the main uncertainties in a limited number 
of unknowns, in order to allow for a simple comparison of the two fluxes. 

We begin by computing the galactic halo emission. This is made by two contributions: one from the 
smooth halo and another by the fraction of the density retained in sub-halos. We closely follow the 
approach of Ref. [l^ , although we retain the angular and energy dependence and use a slightly different 
notation. The differential flux of photons (in units of photons per area per time per steradian per energy) 
from dark matter annihilations with a smooth distribution can be written as 

lUE,^)^ME)uf ds^iHlM^, (1) 
Ji.o.s. 47r 

where 



r{s,i/j) = y Tq + ~ 2rQ scosip, (2) 

■0 is the angle between the direction in the sky and the Galactic Center, tq « 8.0 kpc is the solar distance 
from the Galactic Center, and s the distance from the Sun along the linc-of-sight (l.o.s.). Particle physics 
enters via the term 

f^iE)U^ME)^-^, (3) 

where m-^ is the mass of the WIMP (Weakly Interacting Massive Particle), ((Tannw) is the annihilation 
cross section, and the factor 1/2 enters since we assume that dark matter constitutes its own anti-particle 
(if this is not the case, it should be replaced by 1/4). The function fj{E) is the photon differential energy 
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spectrum per annihilation (with units of E~^). This spectrum per annihilation is largely the result of 
the fragmentation and hadronization of the WIMP's annihilation products. 
A general class of smooth halo distributions can be fitted as 

where pq is the dark matter density at the solar distance from the Galactic Center, and a is a characteristic 
scale radius below which the profile scales as r~'^ . Two of the most well known profiles have been 
proposed by Navarro, Frenk and White (NFW) [J, with 7 = l,a = l,e = 2 and Moore et al. Q with 
7 = 3/2 , a = 1 , e = 3/2. A cored isothermal profile is recovered in the case of7 = 0,a = 2,e = 1 [l5l |. 

These profiles differ primarily in the central region. Since here we are focusing on the Galactic diffuse 
emission rather than that from the Galactic Center, our choice of halo profile is not critical. The uncer- 
tainties which are introduced through the choice of profile (within a factor 2) are negligible for our 
discussion. For definiteness, we shall adopt an NFW profile with /9sm(''o) = 0.3 GeV/cm^ and a = 45 kpc. 

We now turn to the sub-halos contribution to the average fiux. Unfortunately many details enter the 
calculation, including the overall fraction of the dark matter mass which is still in the clumps, the dark 
matter profile of the sub-halos, their mass distribution, etc. However, it has already been shown in the 
past (see e.g. L. Bergstrom et al. in Ref. 01) that their average gamma-ray flux is mainly dependent on 
a single, effective parameter (note that this is not the case for the statistical properties of their angular 
fluctuations). To illustrate this point, let us write the annihilation rate from a single sub-halo, Fci, in 
terms of the density proflle of the clump, Qc\, as 



X ^ 2 gci(x) = nA/go- (5) 



The symbolic evaluation of the integral at the RHS has been written in terms of a mass scale associated 
to the clump, M, and a density parameter, go- In general terms, the flux from the clumpy component 
can be written as 

U{E, i') = ^ ds y dA/ j dp n,i[r(s, V), M; p] T,,{M,p), (6) 

where M is the mass of the sub-halos, p a symbolic variable for all the other unaccounted independent 
parameters, and nc.i[r{s,ip), M; p] is the number density distribution (differential with respect to the M 
and p) of sub-halos as a function of M, p and the distance from the Galactic Center, r. The exact range 
over which the M and p integrals are carried out depends indirectly on the cosmology and the particle 
physics model which is considered. Equation ([5|) motivates one to write 

dM j dpn,,[ris,ij),M,p]T,i{M,p) = go C[?-(s, ^)] /Osm[r(s, V)] , (7) 

where now go is purely a convenient normalization factor, which we flx similarly to Ref. [3l as go = 10~^^ g 
cm~'^. For a constant the RHS of Eq. ([7]) is the mathematical expression of the hypothesis that the 
distribution of dark matter confined in sub-halos traces the smooth distribution of the halo. In this case, 
C, can be thought as measuring the fraction of dark matter mass which is still in the clumps. While 
C ~ const, is approximately true at intermediate galactocentric distances, both semi-analytic arguments 
and simulations suggests that (^{r) Psmir) may tend to a cored distribution at low r, and possibly decline 
at large r as instead of the law expected both from the NFW and the Moore profile fwi. [TtI. fisj . 
Equivalently, one may have approximately that C(r) oc r at large r and low r, while ^ const, at 
intermediate galactocentric distances. For simplicity we shall consider the ( = const, as our baseline 
case, but when relevant we shall comment on the implications of the violation of this Ansatz in the 
following. Also, the distribution of surviving sub-halos is especially uncertain at small galactic radii, 
given the important role that tidal interactions with baryons might play [l^. Since this is still an open 
issue, to limit the impact of these uncertainties, we conservatively discard the predictions of annihilation 
rates in the inner regions r < 3 kpc (|V'| < 20°). The residual effect of interactions with disk stars will 
only affect the overall annihilation rate by a factor of 0(1), which is subleading with respect to other 
unknowns. 

Under the previous assumptions, we can rewrite Eq. ([6|) as 

UE,^)^f,{E)nCgJ d.^^^^- 
Ji.o.s. 47r 
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Combining the smooth and clumpy components, the total galactic signal can be written as 

/gai(-B,^) = + /ci](£;,V^) = ^ME)Il [ dspsmHs,^)] X (p„„[r(s»] +Cgo)- (9) 

Jl.o.s. 

A numerical evaluation shows that even small degrees of substructure (as low a,s ( ^ 10~^) are sufficient to 
cause the sub-halo contribution to dominate over the smooth component everywhere but in the innermost 
region of the Galaxy. It is worth commenting that this conclusion does not only hold for our Galaxy, 
but also for other promising targets for dark matter annihilation signals, like dwarf spheroidal (dSph) 
satellite galaxies of the Milky Way. Although the geometry of the problem is different — in the Galactic 
case we are inside the halo — it can be shown that the dumpiness in dSphs may be responsible for a 
global enhancement with respect to the smooth halo flux of up to a factor ~100 (sec Strigari et al. in 
Ref. [l|). For comparison, a value of C — 0.03 would cause a similar enhancement for the sub-halo to 
halo emission of our Galaxy in the direction perpendicular to the Galactic Plane. Notice that, being qq a 
convenient normalization constant chosen by hand, we have shown that within simplifying but reasonable 
hypotheses the clumpy contribution only depends on the single effective parameter 

On the other hand, the extragalactic contribution to the diffuse gamma ray spectrum is given by [lol | 

i..iE) = ^nf d. j^^^-^ MEH + .)] e-^(-^^) , (10) 

where the Hubble function, H{z), is related to the present Hubble expansion rate, i/o, through the 
matter {^m) and the cosmological constant energy density (^a), H{z) = Ho\/^m{^ + z)^ + ^A- The 
factor e~'^^-^'^^ accounts for the finite optical depth, r, of the universe to high energy gamma rays due to 
scattering with the extragalactic background light. The dumpiness of the dark matter is usually taken 
into account by writing 

^2 /.,\ _ / Pdm,0 \ 2/11 „\6a2/^\ _ n2 I "^^0 



piJz) ^ K~ pUi + zfA\z) = ^ (1 + zf^\zl (11) 

where the so-called enhancement factor, A^(z), can be parameterized approximately as (see, for example, 
Ref. [13]) 

A^(.) = (12) 

^' (1 + Z)3' ^ ' 

where A^(0), at least for z < 10, has only a weak dependence from z. This is confirmed in more detailed 
calculations, as the one in UUio et al. in Ref. [l|. Their Fig. 5 clearly shows that (1 -I- z)^ A^(z)/if(z) cx 
1/H{z), apart for the range z < 1 where a residual dependence of of 0(1) may be present, depending 
on the model considered. Since this is a sub-leading uncertainty for our level of approximation, we 
will neglect any z dependence in A^(0). According to semi-analytical estimates, A^(0) can have values 
ranging from lO"' to 10^ AU together. 

The relative weight of the Galactic and extragalactic contributions has only a weak dependence on the 
spectral shape, In the following, we shall consider the integral quantity, Fi{ip) = H~^ / dEIi{tp). We 
also adopt, for definiteness, = 100 GcV whenever not specified otherwise. 

In Fig. [U we show in the ( — A^(0) plane the three contours i^gaiZ-FLx = 0.1, 1, 10 in the direction 
orthogonal to the Galactic plane (-0 = 7r/2), where the astrophysical contribution to the diffuse flux is 
minimized. For the -Fgai/fcx = 1 case, we also show results for two other values of ip (7i'/4 and 37r/4). 
In the lower-right region of the plot, the Galactic radiation dominates the dark matter annihilation 
contribution to the diffuse fiux (the absolute flux depends on the factor H). To the right of the purple 
vertical line, the contribution from sub-halos dominates over the smooth halo term at -0 = 7r/2. 

Although very rough, the simple model sketched above is useful for understanding qualitatively some 
features, as illustrated by the following example. In Ref. ^l|, the claim is made that a dominant fraction 
of the high energy diffuse gamma ray spectrum measured by EGRET [l^, US] is consistent with being due 
to ^ 500 GeV annihilating neutralinos, without violating existing constraints. Ref. [l^l, in contrast, 
concludes that it is hardly conceivable that dark matter annihilation could be a main constituent of the 
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FIG. 1: The contours -Fgai/Fcx = 0.1, 1, 10 in the direction orthogonal to the Galactic plane {tp = 7r/2) in the 
^ — A^(0) plane, as well as the contour -Fgai/Fcx = 1 for three values of -0. In the lower-right part of the plot, the 
Galactic radiation dominates the dark matter contribution to the diffuse flux. To the right of the purple vertical 
line, the contribution from sub-halos dominates over the smooth halo term at ■0 = 7r/2. 

extragalactic background without exceeding the observed gamma ray flux from the Galactic Center, at 
least if the density profile of the Milky Way is not very different from that found in other galaxies. We note 
here that the two statements arc equally valid within the presently allowed par ameter space, although 
they assume very different priors. The enhancement factor considered in Ref. [2l| is around (0) ^ 10'' 
(where the extragalactic flux is likely to dominate for any reasonable value of C), while in Ref. [2J], values 
of A^(0) ~ lO''^ — 10^ are assumed (for which the galactic flux dominates over the extragalactic one for 
any reasonable value of (). Although the sub-halo profiles assumed are the same, in [2J] it was noted 
that such a large discrepancy may come from uncertainties concerning the concentration parameter and 
the presence/absence of substructures. Whatever the cause, note that already a boost of 10^ — lO'^ in the 
enhancement factor would bring the expectation of Ref. [2J| for the extragalactic contribution within an 
order of magnitude of the EGRET flux, similar to the flndings of Ref. [2l|. Although, as correctly noted 
in Ref. the parameters that we call here C and A^(0) are to some extent correlated, it is plausible 
that a "large" ( value might not affect significantly the Galactic Center signal, but only the Galactic 
diffuse signal, due to the effects of tidal disruption, for example. Namely, the Galactic Center flux, which 
has the highest uncertainty due to the unknown extrapolation of the smooth halo profile, has probably 
no significant contribution from the clumpy fraction of the halo. 

III. DISTINCTIVE PATTERNS OF A DOMINANTLY EXTRAGALACTIC 

CONTRIBUTION 

In this section, we discuss the most prominent features in the angular distribution of the diffuse gamma 
rays if the extragalactic contribution dominates. Although the extragalactic dark matter fiux should 
carry interesting statistical information (see [1^, H^; also, see [13, [13 for a discussion of astrophysical 
backgrounds) we shall concentrate our analysis here on deterministic features, especially at large scales. 

A. Cosmological Compton-Getting Effect 

An observer in motion with velocity, u, relative to the coordinate system in which the distribution 
of gamma rays is isotropic will measure an anisotropic fiux. If gamma ray sources are, on average, at 
rest with respect to the cosmological frame, the magnitude and direction of the velocity, u, of the solar 
system can be deduced from the detection of the dipole anisotropy of the CMB, u = 369 ± 2 km/s in 
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the direction {l,b) = (263.86°, 48.24°) [2^. Since m = |u| ^ 1, the anisotropy is dominated by the lowest 
moment, i.e. its dipole moment. The Galactic analogue of this effect was first noticed as a diagnostic 
tool for cosmic rays by Compton and Getting [soj . 

The amplitude of this anisotropy can be derived from the Lorentz invariance of the the phase space 
distribution function, /, in the frame of the observer and of the emitters (see, for example, Ref. [3l|). A 
first order expansion allows one to deduce the shape of the differential intensity. If we denote with E the 
energy in the emission (rest) frame, and with E' the energy as measured in the observer (moving) frame, 
the intensity in the moving frame can be written as 

I'iE,n)^IiE) 

where in the above formula n is the generic direction and the identification of the energy in the two 
frames E = E' has been made consistently with the first order result. The function, I{E), is the number 
of particles per unit solid angle and unit energy that pass per unit of time through an area perpendicular 
to the direction of observation, and it is related to the phase space density, /, by I{E) ~ E"^ f{E). Thus 
a dipole anisotropy is expected with an amplitude given by 

-^max -^min f d In / 

/max + /mill \ dlvL E 

which is independent of energy as long as the energy spectrum does not change, i.e. when it can be 
approximated by a single power-law. Taking into account the observed spectrum of diffuse gamma rays 
I{E) oc E~'^-^ [ijl, one infers ^ = (2 + 2.1) m ~ 0.5%. Note that the Earth's motion with respect to the 
Sun induces a subleading (8%) modulation in the vector u. 

Unfortunately, this is not characteristic of a dark matter signal, but would be shared by any cosmic 
distribution of sources. Yet, its detection (or limits on its amplitude) would allow one to constrain the 
fraction of the diffuse gamma ray background emitted by cosmological sources [3l| , which would translate 
into a conservative constraint on the product, xlix A^(0). 
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B. Intrinsic Anisotropies 

The distribution of the dark matter in the present day universe is highly structured. As a result, intrinsic 
anisotropies should be present in the diffuse gamma ray background from dark matter annihilations. 
Recentl y, i t has been proposed that one could use the peculiar anisotropy at small-scales to probe dark 
matter [25|. Provided that dark matter contributes in a relevant way to the diffuse flux, this signature 
has promising chances to be detected by GLAST, the Gamma Ray Large Area Space Telescope to be 
launched at the end of this year by NASA. This interesting conclusion was recently extended in Ref. [3^ . 
which notes that most of the anisotropy pattern at relatively large angular scales is due to nearby large 
scale structures. This is especially true at energies higher than 0(0.1-1) TeV, when gammas start to 
be absorbed through pair production on the extragalactic background light, and the observable horizon 
shrinks considerably. This means that one can use the whole information of a real large-scale structure 
catalog, going far beyond mere statistical predictions. How strong the correlation between gamma- 
ray sources and the large scale matter potential is depends on the bias of the sources with respect to 
the catalogue tracers. In general, it is present both for astrophysical and dark matter models for the 
production of gamma rays, although the correlation with overdensity should be stronger for dark matter 
annihilation radiation. For high energy photons, pair production energy loss is sufficiently effective such 
that the anisotropies depend almost entirely on the structure of dark matter in the local universe (within 
several hundred Mpc or so). In this case, cosmological uncertainties, the energy spectrum of the dark 
matter signal (in dark matter models), the redshift evolution of the sources, etc., each play subleading 
roles. Additionally, while the "local" anisotropic flux is almost unchanged, at increasingly high energies, 
the far isotropic component is cut away more and more. The relative anisotropy is higher the greater 
the energy cut. Together with the many modes available {21 + 1 for each multipolc coefficient, C/) and 
the absence of the cosmic variance limitation, this compensates, to some extent, for the lower statistics 
available in the E > 0.1 TeV range. Of course, the strength of this technique depends critically on the 
mass of the WIMP being considered. See Ref. [s^ for further details. 

Additionally, the shrinking of the observable horizon for gamma rays with energies higher than 0(0.1-1) 
TeV has important consequences for the spectral shape of a dark matter signal if m-^ is large enough. In 
the left frame of Fig. [5] we plot the spectral shape of gamma rays produced in dark matter annihilations 
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FIG. 2: Left: The shape of the gamma ray spectrum from dark matter annihilations from four redshifts (0, 0.1, 
0.3 and 1.0), including the effects of redshift energy losses and absorption with the cosmic infrared background. 
In the right frame, we compare the spectral shape of gamma rays from local dark matter annihilations (solid) to 
that from sources integrated out to a redshift of 1, 2 and 10 (from bottom-to-top, dashed). We have assumed that 
A^(0) is roughly constant over this range of redshifts. In each frame, a 500 GeV WIMP annihilating to W^W~ 
was considered. The cosmic infrared background spectrum found in Ref. [s^ was adopted. 



at four redshifts, z =0, 0.1, 0.3 and 1.0. Here, we have considered a 500 GeV WIMP annihilating to 
W~^W~ . In the right frame, we show the shape of the gamma ray spectrum from dark matter annihilations 
integrated out to a redshift of z —1, 2 and 10 (from bottom-to-top, dashed), compared to the harder 
spectral shape emitted at redshift zero, e.g. from a local source (solid). Here, we have assumed that 
A^(0) is constant in this range of redshifts. Also, wc have matched the integrated spectra to the local one 
at Ery = TO^, in order to show clearly how the redshift distribution softens the continuum spectrum. In 
each frame, we have adopted the model of the cosmic infrared background spectrum found in Ref. (33l |. 
The optical depth of the universe to high energy gamma-rays is calculated by simply integrating over the 
cross section for pair production and the spectrum of the infrared background. 



IV. DISTINCTIVE PATTERNS OF A DOMINANTLY GALACTIC CONTRIBUTION 



It is well known that the angular distribution of gamma radiation from dark matter annihilation in the 
Galactic halo has potentially interesting signatures (see, for example, the seminal paper, Ref. 3J| or the 
discussion based on the recent simulations reported in J. Diemand et al. of Ref. Q). In the following, 
we summarize the basic properties expected, once again emphasizing their deterministic features. 



A. Offset Position of the Sun 



In the limit of exact spherical symmetry of the dark matter halo, an observer placed in the center of 
the distribution would observe an isotropic annihilation signal. This follows trivially from Eq. ^ in the 
limit r0 — > 0. However, the Sun is offset with respect to the center of the Galactic halo. This induces 
a peculiar angular dependence, with a maximum toward the inner Galaxy and a minimum toward the 
antigalactic Center. We show this in Fig. [3] for three values of C. In the limit of spherically symmetric 
halo, this is independent of the azimuthal angle, 0, around the Sun-Galactic Center line. This signature 
is fairly robust in two limits: 

A) When ( < 0.0002 (and A^(0) < 10^) and the angular distribution is dominated by the smooth halo 
(bottom curve in Fig. [3]). For example, the ratio of the flux at V' = 45° to the flux a.t ip = 135° is 
about 5. 

B) When C > 0.0005 (and A^(0) < few x 10^) and the angular distribution is dominated by the clumpy 
halo (top curve in Fig [3]). In the limit of a constant the ratio of the flux a.t tp = 45° to the flux 
at V = 135° is of about 2.2. 
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FIG. 3: The angular profile of the Galactic emission of gamma's from dark matter annihilation for three values 
of the dumpiness parameter, ^. The curves have been normalized to 1 at = 20°. The region i/j < 20° is not 
shown, since it is affected by large uncertainties (see text for discussion). 



Note that the two models considered in Ref. [35| correspond to these two extreme situations. It is 



however worth to point out that the quantitative conclusions in the second case (C > 0.0005) are not 
robust if one drops the assumption of a quasi-constant C,{r). In particular, a decreasing (^(r) within 
the solar circle would partially suppress the anisotropy, and paradoxically one may end with an almost 
isotropic flux despite the offset position of the Sun. A quantitative estimate of such effect is unfortunately 
model-dependent, in particular on the radial dependence of the tidal-stripping effects of dark matter sub- 
halos, including the non-negligible effects of the baryonic material in the inner solar system. For details 
on this (still) open issue, we address to Ref. [Toj . 



B. Compton-Getting Effect (Proper Motion of the Sun in the Halo) 



A second signature, which to the best of our knowledge is discussed in detail here for the first time in 
relation to dark matter, is purely kinematic. Unlike the disk of our Galaxy, which is supported against 
radial collapse by its angular momentum, the dark halo is supported by random velocities which serves as 
a coUisionless pressure. It is expected that the velocities of dark matter particles were isotropized at the 
time of the formation of the Galaxy in the so-called process of "violent relaxation" [s^ , and should retain 
this distribution as long as no relevant interaction with the collisional baryonic gas intervenes (for a more 
extended discussion, see Ref. [s^]). Although the halo may have some rotation, we know for sure that it 
is not supported by angular momentum, otherwise it would be flattened similarly to the disk. For our 
Galaxy, observational constraints imply that the halo is almost spherical, with differences between the 
axes of the best- fit spheroid not larger than ~ 20% [s^ [s^ . If we neglect an intrinsic angular momentum 
of the halo, our motion with respect to the halo is fully due to the Galactic disk rotation around the 
center of the Galaxy with a velocity of about 220 km/s at the solar galactocentric distance (see e.g. 
Ref. [3§] and references therein). Analogously to the previously discussed extragalactic case, this would 
be manifest as a halo Doppler effect. In this Galactic case, the dipole points in the direction of motion of 
the rotation of the Galaxy, and with an amplitude of about 0.3% [see Eq. ([T5|)]. The amplitudes increase 
when the spectrum softens, and may be very large (although in a narrow energy range) close to an abrupt 
cutoff in the energy spectrum, as the one for Ey ~ m^. Note that in the case of a Galactic astrophysical 
(non-dark matter) origin of the diffuse gamma ray flux, this effect is not present, since the disk is on 
average co-moving with the Sun, with typical proper velocities much smaller that the bulk one. This 
signature is, therefore, specific to dark matter, unlike its extragalactic counterpart described in Sec. lIII A1 
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C. Asphericity of the Halo 

Thus far, we have imphcitly assumed that the dark matter distribution in the halo has a spherical 
symmetry. Simulations of the halo structure, however, show some departure from the spherical symmetry, 
and typical dark matter halos are actually spheroidals with a moderate eccentricity. Also, the level of 
ellipticity may be radially varying, with larger ellipticities towards the inner Galaxy, see e.g. [i^. Note 
however that there are hints that the inclusion of baryonic dissipation may produce significantly rounder 
halos than those formed in equivalent dissipationless simulations especially in the inner halos [41| . 

The impact of this effect on the annihilation flux was discussed in detail in Ref. [s^l, and we shall 
not repeat here that analysis. Qualitatively, the anisotropy in the emission may be comparable to the 
effect of the offset position of the Sun. However the symmetry axes of the halo are unknown. Thus 
the magnitude, the exact angular dependence and the direction of this anisotropy are unknown, making 
identification of this effect challenging. 

D. Proper Anisotropy Due to dumpiness 

As in the case of extragalactic emission, the typical parameters of the clumps are reflected as peculiar 
angular signatures of the annihilation signal in the sky. In contrast to the signatures considered above, 
this is most important at small scales, unless a few clumps happen to be very near us (and in the case 
that they are resolved, it may also be possible to detect their proper motion (12). While a highly 
clumpy and steady emission may be a spectacular signature of dark matter annihilation, it is difficult 
to obtain model-independent predictions on the expected angular power spectrum, since the physical 
processes determining the result are highly non-linear and entangled (such as mergers, tidal disruptions, 
etc.). Furthermore, the observer-related "galactic variance" — the dependence of the signal features on the 
actual distribution of clumps around our position, which is unknown — may greatly change the expected 
signal. Some attention has previously been given to this kind of signature within the context of superheavy 
dark matter decays connected with the ultra- high energy cosmic rays j43j . 

V. PROSPECTS FOR GLAST 

Regarding the detectability of dark matter annihilations through the various angular signatures de- 
scribed in this paper, we shall limit our discussion to the GLAST satellite detector [3]. Although the 
curre nt g eneration of atmospheric Cerenkov telescopes, including H.E.S.S. [1^, MAGIC [i^ and VERI- 
TAS [43], are highly sensitive to very high energy gamma rays, possibly including dark matter radiation, 
they are best suited for studying point sources rather than diffuse emission. 

The flux of gamma rays due to dark matter annihilation is at present unknown. To the best of our 
knowledge, it may account for a signiflcant or even dominant fraction of the diffuse gamma ray back- 
ground, especially at the highest measured energies [2T|. A residual, isotropic radiation usually interpreted 
as extragalactic gamma ray background has been measured by GLAST's predecessor, EGRET j22||. The 
intensity of this spectrum can be flt by [l^ (see also Ref. (23j l 

/ fj \ -2.10±0.03 

/cosmic(-B^) = (7.32 ±0.34) X 10-6 cm-\-hT-'GcY-\ (16) 

over an energy range from ~ 10 MeV to ^ 100 GeV. The number of events collected above an 
energy, E-y, can be written as 

poo 

N^ = t- rifov • / dE A,s{E)I,osmi,{E) , (17) 

where fifov is the solid angle of the field-of-view, Aos{E) is the effective collecting area of the instrument 
(averaged over its fleld-of-view) , and t is the time observed. The LAT detector on board GLAST features 
Aeff (-B) ~ lO'' cm^ iorE>l GeV (and Aes{E) > 4000 cm^ for E > 0.1 GeV), 17fov 2.4 sr, and excellent 
hadronic background rejection capabilities |44l ]. Even when considering realistic fractions of data to be 
rejected from the use of mask cuts or a non-ideal duty-cycle, at least 0(10^) diffuse photons per year 
should be measurable in the GeV range, given the background flux measured by EGRET. Considering 
up to a decade of operation, it is clear that as long as dark matter is responsible for a significant fraction 
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FIG. 4: Contours in the — A'^(O) plane (see Sec. |TTJ illustrating the fraction of the diffuse flux observed by 
EGRET (above 1 GeV) which is the product of dark matter annihilations. Above and to the right of the solid 
line, more than 10% of this flux is produced by dark matter. Above and to the right of the dashed line, dark matter 
generates 1% or more of this flux. Here, we have used a 100 GeV WIMP with a 3 x 10~^^ cm'^/s annihilation 
cross section. The Galactic flux being considered here is at a direction perpendicular to the Galactic plane. 



of the flux of Eq. ([16]), say at least 10-15%, GLAST should be sensitive to anisotropies down to the few 
permil level. Modulo the issue of foreground removal, this precision should be sufficient to detect the 
tiniest of the signatures considered in this paper. For illustrative purposes, in Fig.|3]we show the regions 
in the parameter space where the DM flux amounts to more than 10% or 1% of the EGRET flux given 
by Eq. ([16]), assuming typical parameters for the DM candidate, and fluxes integrated above 1 GeV. 
In the region above the solid line the collected statistics would allow one to detect all the signatures 
discussed in this article. Most of these signatures should be detectable in the region above the dashed 
line. Interestingly, in the lower left corner of the parameter space where the DM flux is subdominant, the 
signatures of the smooth halo signal should be prominent. Since the offset position of the Sun implies 
anisotropies of 0(100%), even in this range the prospects for the detection of DM signatures are more 
promising than what one would naively deduce from Fig. [U 

Of course, a major obstacle in detecting such features are astrophysical foregrounds, which can make 
the identification of the above signatures difficult (the same is true, of course, when considering spectral 
energy features). The very subtraction of the Galactic astrophysical foreground is a highly non-trivial 
issue, see e.g. [23|. For the purpose of detection of the diffuse dark matter flux, apart for applying low- 
latitude cuts (the astrophysical contamination is higher at low galactic latitudes), it may be appropriate 
to fit simultaneously for the sum of Galactic astrophysical ffux (as provided by numerical routines like 
GALPROP [4§|) plus an additional dark matter component. In addition, there are fortunately a number 
of qualitative properties of dark matter signals which may help in distinguishing them from astro phy sical 
emissions. The strategy to reveal peculiar extragalactic features has been discussed elsewhere (26l.l32l| and 
we shall not repeat it here. We want only to stress that, qualitatively, after bright sources are removed, a 
significant cross-correlation with large scale structure catalogs or the CMB dipole direction would provide 
important diagnostics and a tool for background rejection. 

In case of dominance of the Galactic dark matter halo emission, it is important to note that: (i) The 
offset position of the Sun produces a ■(/'-dependent asymmetry, (ii) The Doppler effect due to our motion 
in the halo produces a dipolar asymmetry in the azimuthal angle, (f>, which for a fixed '0 identifies the 
direction around the Sun-Galactic Center axis (each corona should share the small dipole in the direction 
of the motion of the disk in the halo). Note also that, although small, its angular shape and direction are 
known a priori, which facilitates a search for it. (iii) The signature of the asphericity of the halo depends 
on the unknown shape of the halo. On general grounds, one should expect a modulation of the emission 
in the angle, ip, (as does the one treated in Sec. IIV Ap . as well as a modulation in the azimuthal angle, 
(as docs the one treated in Sec. IIVBI which, however, should be smaller). Barring fine-tuned directions 
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of the halo axes of symmetry, a more specific signature is an asymmetry pattern between the 4> 
—ip, (p directions (reflection with respect to the Galactic plane), which cannot be mimicked by any of the 
two previously considered effects. Note that, although the relative magnitude of the anisotropy measured 
may change (whenever dark matter does not constitute a constant fraction to the diffuse gamma signal), 
the relative weights of the different Galactic asymmetries discussed here remain constant with energy. 
Since the diffuse radiation will be observed in many energy bands, one should keep in mind that the 
shape of this pattern must stay constant at each energy if it is to be attributed to dark matter. Also, if 
suspected dark matter spectral features appear, they should correlate with the amplitude of the angular 
asymmetries discussed in this paper. Combining angular and energy spectral information should surely 
strengthen our ability to distinguish dark matter annihilation radiation from other types of sources. 
Indeed, there is no a priori reason to expect that an astrophysical foreground affecting cither the energy 
spectrum or the angular spectrum should appear in the other as well. 

VI. CONCLUSIONS 

With the next generation satellite-based gamma ray telescope, GLAST, our measurements of the diffuse 
gamma ray spectrum will improve dramatically. Not only will the number of events observed increase, 
but also the spectrum will be studied up to considerably higher energies. Furthermore, much of the 
astrophysical contributions to the diffuse flux observed by EGRET will likely be resolved as point sources 
result of GLAST's superior angular resolution. 

In this article, we have studied the signatures and distinctive features that would be possessed by 
dark matter annihilation radiation from cither a dominantly Galactic or cxtragalactic population of dark 
matter. For an cxtragalactic contribution, we discussed the effect of observing dark matter radiation in 
a moving frame of reference (relative to the dark matter distribution), known as the Compton-Getting 
effect. We also discuss the ability to tie the angular distribution of annihilation radiation to the known 
structure of our (cosmologically speaking) local universe. This is especially interesting for gamma rays 
more energetic than ~100 GeV, which can be absorbed by the cosmic infrared background radiation, 
effectively reducing the horizon for such particles. 

If the diffuse dark matter annihilation radiation is instead dominated by the Galactic population, this 
flux will contain distinctive features resulting from the offset location and proper motion of the Solar 
System relative to the Galactic halo. Furthermore, any asphericity in the halo will have effects of the 
angular distribution of annihilation radiation. As in the extragalactic case, we also briefly discuss the 
anisotropics which are likely to result from inhomogeneities and substructure in the Galactic dark matter 
distribution. 

GLAST's predecessor, EGRET, observed the presence of a diffuse gamma ray spectrum. Based on 
these observations, GLAST should see O(IO^) diffuse photons per year above the GeV ^. If a signiflcant 
fraction of this flux is the product of dark matter, many of the angular features described here could 
potentially be identified. For example, if 10-15% of the diffuse fiux observed by EGRET is dark matter 
annihilation radiation, in a decade of operations GLAST should be sensitive to anisotropics down to the 
few permil level, comparable or often better than the level required to study the angular signatures we 
have discussed in this paper. 



Acknowledgments 

This work has been supported by the US Department of Energy and by NASA grant NAG5-10842. 

APPENDIX A: FLUX ENHANCEMENT DUE TO CLUMPINESS 

Consider a cubic volume of size L at a distance D » L from us, such that the solid angle subtended by 
its surface is O^, where Q = L/D. Let us assume that it is filled with: (i) a homogeneous distribution of 



At high latitudes and for energies E > GeV, the diffuse EGRET fiux is comparable or larger than the Galactic flux 
reported in l49t . as it can be easily checked with the fit provided by L. Bergstrom et al. in Ref [l]]. 
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dark matter with density p; (ii) a cubic dark matter clump of size I <^ L and density pc, having the same 
mass {i.e., such that pc = pL^)- We denote with / and Ic the differential photon fluxes (per area, per 
time, per solid angle) collected in the two cases, and with J = IQ^ and = Ic the flux per area per 
time, where 9 = l/D. If the angular resolution of our instrument is larger than 0, then the quantities J 
and Jc are the relevant observables. 

Let us assume now that / oc /j ^ ^ p'^dx, where k = 1 corresponds to a decaying dark matter scenario 
and K = 2 to an annihilating dark matter one. One immediately sees that 

Ic<^P:i^P^(^j^'^l = P^(^jy^ \, (Al) 

and 



3k-1 



(A2) 



from which it follows 



Equation (|A3p implies that when k = 1 (the decaying dark matter scenario) , only the angular information 
is sensitive to the dumpiness of the distribution. In other words, the total photon flux from decaying dark 
matter can be calculated (for a fixed particle physics scenario) once the profile of the average dark matter 
density is specified. This implies for example that the Galactic halo signal dominates over the cosmological 
one. For the ti ~ 2 (annihilating dark matter) scenario, however, the picture changes considerably. In 
particular, the dumpiness of the dark matter distribution strongly affects both the angular distribution 
of the emission and the overall flux. 
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